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Abstract. A model, plane symmetric, 3-D potential, which 
preserves some features of galactic problems, is used in or- 
der to examine the phase space structure through the study of 
the properties of orbits crossing perpendicularly the plane of 
symmetry. It is found that the lines formed by periodic orbits, 
belonging to Farey sequences, are not smooth neither continu- 
ous. Instead they are deformed and broken in regions charac- 
terised by high Lyapunov Characteristic Numbers (LCN's). It 
is suggested that these lines are an incomplete form of a frac- 
tal repeller, as discussed by Gaspard and Baras (1995), and 
are thus closely associated to the "quasi-barriers" discussed by 
Varvoglis et al. (1997). There are numerical indications that the 
contour lines of constant LCN's possess fractal properties. Fi- 
nally it is shown numerically that some of the periodic orbits 
-members of the lines- belong to true irregular families. It is 
argued that the fractal properties of the phase space should af- 
fect the transport of trajectories in phase or action space and, 
therefore, play a certain role in the chaotic motion of stars in 
more realistic galactic potentials. 



Key words: galaxies : kinematics and dynamics - chaos - stel- 
lar dynamics 



1. Introduction - Motivation 

One of the most interesting open questions in non-linear dy- 
namics is the nature and evolution of transport in the chaotic 
phase space regions of perturbed integrable Hamiltonian sys- 
tems. Most of the work in this field has been done for 2-D 
systems or for the standard map, following two different ap- 
proaches. In the first approach, which may be designated as 
microscopic, one studies the homoclinic and heteroclinic tan- 
gle of stable and unstable manifolds inside the chaotic region. 
This gives a "complete" qualitative picture of transport in phase 
space, since a chaotic trajectory follows the unstable manifolds 
of the unstable periodic orbits, jumping from one to another 
near the heteroclinic points. On the other hand, invariant tori 
surrounding stable periodic orbits act as "barriers" in the trans- 
port process, modifying at the same time the topological struc- 



Send offprint requests to: H. Varvoglis 



ture of phase space. However, it is not presently well under- 
stood how the results from this kind of study can be used in 
a quantitative description of the evolution of an ensemble of 
trajectories. 

In the second approach, which may be designated as macro- 
scopic, one assumes a priori that transport is a pure random 
phenomenon (a Markovian process). Then the transport of an 
ensemble of trajectories may be considered as "classical diffu- 
sion" in action space. Unfortunately, it is well known that trans- 
port in Hamiltonian systems cannot be considered as a pure 
Markovian process, due to the problem of stickiness. More- 
over the process may follow Levy rather than classical statistics 
(Shlesinger et al., 1993). In this case the macroscopic approach 
still works, provided that the process is considered as fractal 
diffusion in normal space (Zaslavsky, 1994). In this way it be- 
comes evident that Levy-like statistics are closely related to the 
self-similarity of the phase space. This formalism results in a 
differential equation with fractional partial derivatives, whose 
solution is formidably difficult, even in the simple case of 2-D 
systems. 

There are rather few attempts to study transport in Hamil- 
tonian systems with more than two degrees of freedom, mainly 
because the above mentioned approaches for 2-D systems can- 
not be directly generalized. It is therefore important to note that 
there exists yet a third approach for the description of trans- 
port, which may be implemented in a straightforward way in 
systems with more than two degrees of freedom, and this is 
normal diffusion in fractal space (e.g. see West and Deering, 
1994, and references therein). The fractality of phase space has 
been already demonstrated for 2-D systems and has been as- 
sociated with the self similarity of the hierarchical structure of 
island families on a surface of section (e.g. see Zaslavsky, 1994, 
BenKaddaetal., 1997). 

In the case of more than two degrees of freedom the sit- 
uation is more complicated, since the topological structure of 
phase space in any region depends on the number of local in- 
tegrals of motion. Recently Varvoglis et al. (1997) have found 
indications of phase space fractality in the same model 3-D 
Hamiltonian system studied in the present work, but they were 
unable, due to the particular method of study they had selected, 
to identify actual fractal or multifractal sets in phase space. 
They have conjectured that the fractality is due to the presence 
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of quasi-barriers (Varvoglis and Anastasiadis, 1996), i.e. ge- 
ometrical objects of lower than full dimensionality (depending 
on the number of local integrals of motion), such as periodic or- 
bits, invariant tori around them and cantori (Aubry, 1978; Per- 
cival, 1979). 

From the above discussion it is obvious that the importance 
of the fractality of phase space in the diffusion process lies in 
the nature and in the rate of the diffusion: In a simply connected 
space, the diffusion is classical and follows Fick's law (trans- 
port proportional to In a fractal space the diffusion is 
non-classical (presumably a Levy process) and, in the case of 
a Hamiltonian system, usually it has a lower than the Fickian 
rate. 

On the other hand, the distribution of periodic orbits per- 
pendicular to the x — y plane in the present 3-D model po- 
tential has been studied by Barbanis and Contopoulos (1995) 
(henceforth referred to as BC) and Barbanis (1996). It was 
found that it is not random: the perpendicular crossings (p.c.'s) 
of orbits with multiplicities forming Farey sequences with the 
x — y plane are arranged along lines, named in the present pa- 
per basic-lines (BL) and Farey tree lines (FTL). These lines are 
not continuous but they possess gaps; two of them have already 
been mentioned in BC. Since the presence of Farey sequences 
implies a form of self-similarity, it is possible that the BL's and 
FTL's may be related to the above-mentioned quasi-barriers 
and that the gaps may be related to any fractal properties of 
phase space. 

In the present paper we use the model potential studied by 
BC in order to test the conjecture put forward by Varvoglis et al. 
(1997) about the origin of the fractality of phase space. Further- 
more, through this, we attempt to assess the relation between 
the topological structure of the phase space and transport, on 
the one hand, and the system of BL's and FTL's, on the other. 
Although this potential can model a true galaxy only locally, its 
study is expected, nevertheless, to contribute to the understand- 
ing of galactic evolution problems, based on the fact that evo- 
lution is intimately connected to transport in phase and config- 
uration space. Transport in a chaotic region of a perturbed inte- 
grable dynamical system, in turn, may be viewed as (stochastic) 
diffusion, which depends on the value of the Lyapunov Char- 
acteristic Number as well as on the topological structure of the 
phase space. We believe that similar behaviour would charac- 
terise any 3-D perturbed integrable dynamical system, as are 
most of the realistic model galactic potentials. 

The purpose of this paper is threefold: 

(a) to study in finer detail the distribution of the periodic 
orbits along the BL's and FTL's. 

(b) to examine whether these lines are connected, in any 
way, to transport in phase space and are, thus, related to the 
"quasi-barriers" and 

to confirm the existence of true irregular families of peri- 
odic orbits in a 3-D Hamiltonian system. 

This paper is organized as follows. In the next section we 
present, for reasons of completeness, some remarks on previ- 
ous results. In Section 3 we show that the LCN's may "map" 
the fractality of phase space. In Section 4 we present our results 



on the interconnection and discontinuities of BL's and FTL's. 
In Section 5 we show numerically the existence of irregular 
families of periodic orbits. Finally in Section 6 we summarise 
and discuss our results. 

2. Remarks on previous results 

2.1. Basic lines 

The systematic exploration of the phase space structure of 3- 
D Hamiltonian systems was initiated by Magnenat (1982) and 
later by Contopoulos & Barbanis (1989), BC and Barbanis 
(1996), through the study of periodic orbits of unit mass test 
particles in a 3-D model potential, corresponding to the Hamil- 
tonian 

h= Upl+pl+pD + hi^ + By' + cz 2 ) (i) 

— exz 2 — r\yz 2 = h 

with a particular set of parameters (A = 0.9, B = 0.4, C = 0.225, 
e = 0.560, 7] = 0.20 and h = 0.00765). The specific values 
for the parameters A, B, C were selected by Magnenat (1982) 
because the 2-D system x — z for y/ A/C = 2 is topologically 
equivalent to the Inner Lindblad Resonance and the 2-D system 
y — z for \J B jC = 4/3 is a well studied 2-D dynamical system 
with extended chaotic regions. For this set of parameters it was 
found in BC that the distribution of the p.c.'s of the periodic 
orbits is not accidental: the crossings are arranged along partic- 
ular lines on the x — y plane (where x — A 1 / 2 x,y = B 1 l 2 y), as 
it is evident in Fig. [j] (taken from Fig. 5 of BC). Each orbit has 
either one p. a, marked by a (x), or two p.c.'s, marked by (+). In 
this paper, for clarity reasons, we use a point (•) instead of (x). 
Each p.c. is designated by a number (or a number and a letter), 
representing the multiplicity of the orbit; different orbits with 
the same multiplicity are differentiated by a prime. The lines 
are formed by the p.c.'s of the orbits, whose multiplicities form 
arithmetic progressions with increment the multiplicity of the 
limit orbit. They are termed in this paper basic if they connect 
orbits of the unperturbed system (e.g. la, lb in Fig. 1) or orbits 
with low multiplicity, m (e.g. m < 5). For example the basic 
sequences 

3(+),5(-), 7(+), 9(0, ...2c(+) 
and 3(+), 5'(-), 7(+), 9'(-), ...2c(+), 

which form the lines C and C', start with an orbit of multiplicity 
3 and increment 2, which is the multiplicity of the limit orbit 
2c. 

2.2. Farey tree lines 

In BC it was recognized that the distribution of p.c.'s on the 
x — y plane presents some self-similar features. To begin with, 
to each basic sequence of orbits corresponds a number of new 
sequences of higher order. Each member of a new sequence is 
a periodic orbit having as multiplicity the sum of the multiplic- 
ities of two consecutive orbits of the generating sequence. For 
example between the first two orbits of the basic sequence C, 
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i.e. 3(+) and 5(0, there is the orbit 8(+). Two second order se- 
quences are formed between 3 and 8 as well as between 8 and 
5, i.e. 

8(+),ll(-),14(+),17(-),20(+), 23(-), ...3(-) 
and (2) 
8(+),13(+),18(+),23(+),28(+),...5(0 

Both new sequences have as first member the orbit 8, but they 
have different limit orbits and increments (i.e. 3 for the first and 
5 for the second, respectively). In the same way we can form se- 
quences of even higher order. We call the orbits of these higher 
order sequences Farey tree orbits and the lines on which they 
lie Farey tree lines, because such sequences are similar to the 
Farey sequences discussed by Niven and Zuckerman (1960). 
For a good review on Farey trees see Efthymiopoulos et al. 
(1997). 

The periodic orbits are used, in the present paper, as an ad- 
ditional tool for probing the topological structure of the phase 
space. In particular, periodic orbits of high multiplicities are 
needed in order to search for self-similar structures in the BL's 
and FTL's for a wide interval of multiplicities (see Avnir et al., 
1998). The importance of the non-basic (Farey) lines lies ex- 
actly on the above line of reasoning, as well as in the fact that 
they might be necessary in order to calculate correctly an ap- 
propriate diffusion coefficient. It should be noted that we have 
restricted our study to the periodic orbits crossing perpendicu- 
lar the x-y plane for simplicity reasons. 

2.3. Irregular orbits 

All periodic orbits of the same multiplicity that can be found 
through the continuous variation of a parameter in the equa- 
tions of motion belong to the same family. The graph giving 
one of the initial conditions of the periodic orbits as a function 
of the parameter is the characteristic of this family. Families, 
which do not exist for values of the perturbing parameter be- 
low a critical value, appear only as pairs, are well known in 
2-D dynamical systems and are called irregular families (Con- 
topoulos 1970; Barbanis 1986). However, their existence has 
not been proven for 3-D systems. Pairs of families that have 
been found in such systems without a direct connection to a 
basic orbit, they were all of small multiplicities and it turns out 
that they have an indirect connection to a basic orbit through 
other families bifurcating from it (Barbanis, 1996). It is desir- 
able to know whether families of periodic orbits in 3-D sys- 
tems (presumably of high multiplicities) do exist without any 
connection to a basic orbit. Their presence might play a certain 
role in the fractality of phase space and, therefore, in trajectory 
transport. 

2.4. Lyapunov Numbers 

Varvoglis et al. (1997), while studying numerically transport 
phenomena in the trajectories of the Hamiltonian (1), have 
found, in an indirect way, that the phase space shows a frac- 
tal structure. In particular they presented numerical evidence 
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Fig. 1. The majority of the perpendicular crossings of the peri- 
odic orbits of the Hamiltonian (1) lie on particular lines on the 
x — y plane inside the boundary circle x 2 + y 2 = 2h. Most 
of the lines are concentrated around the basic orbits la and lb. 
There is also a noticeable spiral structure with focal point the 
orbit la. Four spirals illustrate the area of this structure. Note 
that in this figure (reproduced from BC, Fig. 5) a different nota- 
tion is used than in the rest of the present paper. Each symbol, 
(x) or (+), shows a periodic orbit with one or two crossing 
points, respectively. 

that the function dVj dt has multifractal properties, where V (t) 
denotes essentially the coarse-grained volume of phase space 
visited by a trajectory up to time t. They found, also, that the 
evolution of the function dVj dt is closely related to the evolu- 
tion of the function 



. . 1 d(t) 



(3) 



whose limit, as t — > oo, is defined as the LCN. The function 
x(t) shows plateaus in the time intervals where dV/dt is close 
to zero and steeply increasing transient segments in the time 
intervals where dV/dt takes large positive values. This fact 
shows that a trajectory is confined successively in regions of 
phase space where the LCN converges to a limit and, there- 
fore, the fractal properties of dVj dt should arise from the (ap- 
parently) self-similar distribution of periodic orbits and other 
quasi-barriers in phase space (Zaslavsky, 1994). 

It should be emphasised that, in the present paper, LCN's 
are not calculated in order to estimate "rates of diffusion", since 
LCN's alone cannot describe completely this phenomenon in 
the case where there exist invariant tori of considerable mea- 
sure. We do, however, calculate LCN's in order to "probe" the 
phase space structure, i.e. to delineate, in an independent from 
other methods way, the "topology" of the phase space. This 
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Fig. 2. A segment of the BL between the p.c.'s of the periodic 
orbits 3 and 5 (indicated by a bold line) and the network of 
various FTL's (corresponding to the upper central part of Fig. |], 
line C). The nonexistent part of C between 8 and 14 is indicated 
by a dashed line. The contour lines of constant LCN's A = 
10 -4 are superimposed and the regions of A > 10~ 4 are gray 
shaded 

new method is based on the property that trajectories starting 
in a certain region, bounded by "quasi-barriers", are charac- 
terised by a certain "Local Lyapunov Number". By drawing 
the contour lines of LCN's, one draws, essentially, the geomet- 
ric boundary of the region bounded by the "quasi-barriers". 

3. LCN's and fractal properties 

The "classical" calculation of LCN's requires a continuous in- 
tegration of the corresponding trajectories, until the function 
x(t) has reached a plateau. It is, therefore, obvious that this 
method cannot be used in a "mass production" procedure. For 
this reason we decided to use a different approach, by calculat- 
ing the value of the function x(t) at various integration time in- 
tervals, At, up to At = 3 10 5 . Comparing the results for differ- 
ent At's we have found that they are qualitatively the same (i.e. 
they form level lines with increasing complexity as one goes to 
finer details), provided that At > 10 5 . So, even if \{t) does not 
reach a plateau, its value at the end of the corresponding time 
interval At > 10 5 gives a correct estimate, at least to order 
of magnitude, of the "degree of stochasticity" of the trajectory 
in the phase space region restricted by the quasi-barriers. We 
designate this value as the Fixed Time LCN. In what follows 



we use, for convenience, the notation LCN. Through this kind 
of LCN's we attempt to probe the topological structure on the 
x — y plane as follows: 

We draw a dense mesh on the x — y plane (for < x < 
0.020 and 0.010 < y < 0.045 at intervals Ax = Ay = 0.0005) 
and we use the nodes of this mesh as initial conditions for tra- 
jectories starting perpendicular to this plane. We calculate the 
LCN of each trajectory at the end of various time intervals and 
we draw the contour lines of constant LCN on the x — y plane 
at the LCN value A = 10~ 4 . In Fig. | we show our results 
for the LCN's at the end of a time interval At = 3 10 5 . Fig. § 
contains only a small part of the studied region. The structure 
of the contour lines in the whole region as well as their varia- 
tion with At is still under investigation and will be the topic of 
another publication. 

The emergent picture is very interesting. One can immedi- 
ately see that the studied area may be divided into two regions: 
Most of it is characterized by small values of LCN's (white re- 
gions); inside the white regions one can find small elongated 
regions characterized by high values of LCN's (gray regions). 
It should be noted that the position and elongation of these re- 
gions is closely related to the position and the direction of BLs 
and FTL's (see next Section). This fact may be interpreted in 
the following way. Since a considerable fraction of the periodic 
orbits along these lines are unstable, the exponential divergence 
of trajectories in the area around them is governed mainly by 
the unstable manifolds of the orbits-members of the lines. The 
situation is reminiscent of a fractal repeller as defined in Gas- 
pard and Baras (1995), i.e. a set of countably infinite unstable 
periodic orbits, which has zero Lebesgue measure and a finite 
Hausdorff dimension. In our case the orbits-members of BLs 
and FTL's may be considered as forming a sort of an incom- 
plete (since not all periodic orbits are unstable)/racfa/ repeller. 
In other words the BLs and FTL's are closely associated with 
the "quasi-barriers", which are geometrical objects that sepa- 
rate "loosely" different phase space regions with different Lo- 
cal LCN, as it has been discussed by Varvoglis et al. (1997) and 
Tsiganis et al. (1998). 

Since (I) there is evidence that BLs and FTL's have frac- 
tal properties and (ii) there is a relation between the BLs and 
FTL's, on one hand, and the LCN's on the other, it should be 
interesting to examine whether the contour lines have fractal 
properties as well. To do so we calculate the LCN's in the re- 
gion -0.02 < x < 0.01,0.05 < y < 0.10 using a coarse 
(Ax = Ay = 0.005) and a fine mesh (Ax = Ay = 0.001). 
This region was selected because it was studied, although with 
a much coarser mesh ( Ax — Ay = 0.01) by Contopoulos 
and Barbanis (1989, Fig. 9). We then plot both contour lines at 
A = 10~ 4 , on the same graph, as shown in Fig. |3|. It is easy 
to see that, while the dashed line, corresponding to the coarser 
mesh, seems rather smooth, the continuous one, correspond- 
ing to the finer mesh, shows "tongues" that oscillate about the 
dashed curve, a picture reminiscent of the classical examples 
of fractal curves. Of course this cannot be taken as a proof, 
even numerical, that contour lines are fractal curves, since the 
available data are not exhaustive (see also Avnir et al., 1998). 
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However, we believe that Fig. |J if considered together with 
the fact that LCN's reflect the distribution of the unstable man- 
ifolds of unstable periodic orbits, is a strong indication that the 
contour lines have, indeed, fractal properties. 



4. BL and FTL are not simple 

As it is evident from Fig. [l] and the relevant discussion in BC 
and Barbanis (1996), the BL's formed by the lowest order se- 
quences look continuous and smooth. Because of this it had 
been assumed that the lines appearing in Fig. |l] contain all the 
higher order sequences described in Section 2 as well, so that 
any fractal properties, arising from the Farey tree character of 
the p.c. sequences, are restricted along the lines and not across 
them. However, the existence of two small gaps in a BL, no- 
ticed by BC, gave the first evidence that, somehow, this picture 
is not correct. In order to examine, therefore, in depth the situ- 
ation, we proceed in the calculation of periodic orbits belong- 
ing to several higher order sequences than those appearing in 
Kg. I 

We focus our interest on a small region, between the orbits 
3(+) and 5(-) of the line C (BC, Fig. 8), since it is in this area 
that the first noticeable gap was observed. Between these two 
orbits lies the orbit 8(+). Therefore two second order sequences 
are formed with orbit 8 as first member and increments 3 and 
5 respectively, i.e. the sequences (2). In Fig. || we can see that 
the first sequence presents a large gap between orbits 8 and 
14, where higher order FTL's, branching away from the BL, 
may be observed. Furthermore, as we discuss below, orbit 1 1 
(which lies between orbits 8 and 14, instead of being part of 
the BL) belongs to a higher order FTL. In contrast the sec- 
ond sequence does not show any noticeable gaps. The reason 
for this difference becomes obvious, if one considers the BL's 
and FTL's in association with the degree of stochasticity of the 
phase space where these orbits lie, estimated through the calcu- 
lation of LCN's. In Fig. || one can immediately notice that the 
first sequence crosses a region of high LCN's, while the second 
is confined in a region of low LCN's. 

Let us now examine what happens to the higher order se- 
quences between 11 and 8, on one hand, and 11 and 14 on 
the other, i.e. ll(-), 19(+), 27Q, 35(+), 8(+) and 14(+), 
25(+), 36(+), ll(-). As it is evident from Fig. f| the FTL 
corresponding to the first sequence is torn and split into two 
branches, one beginning from orbit 1 1 going upwards and the 
other beginning from orbit 8 going downwards. Similar is the 
situation with the second sequence. The FTL is torn and split 
into two branches, one going upwards from orbit 14 and the 
other going downwards from orbit 1 1 . The branch going down- 
wards from orbit 11 is connected to the FTL of the first se- 
quence going upwards from orbit 11, resulting in a compos- 
ite S-shaped line (the three leftmost continuous thin lines in 
Fig. ^J). In this way we see that the gap between orbits 8 and 
14 is, in a way, "filled" with FTL's corresponding to higher or- 
der sequences. This is an example of a case where higher order 
FTL's do not lie on the "parent" BL. 



It is interesting to note that the branch of the second FTL 
going upwards from orbit 1 1 is developed inside a narrow strip 
of ordered motion. 

In Fig. H we notice also two "vertical" higher order FTL's 
that emanate from the orbits 57(-) and 17(-), corresponding to 
higher order sequences with increment 40. The first of them, 
the one emanating from orbit 57(-), is vertical to the local gra- 
dient of the contour lines and lies in a region of low LCN's. 
The second one, the one emanating from orbit 17(-), is vertical 
to the local gradient of the contour lines as well but it seems 
to span a high LCN's region. The most probable explanation is 
that the line lies, in fact, in a narrow strip of low LCN's, with a 
width less than the size of the mesh used. This is something that 
should be expected for a geometrical object with fractal prop- 
erties and has been indeed encountered in some other regions 
of the studied area. 

From the above discussion emerges an intuitively appealing 
picture: BL's and FTL's are, in a sense, lines representing some 
ordered features of the dynamical system, which are deformed 
in the neighbourhood of the chaotic seas and tend to run per- 
pendicular to the local gradient of the contour lines of constant 
LCN's. Therefore we understand that the structure of the BL's 
and FTL's on the x — y plane is considerably more complicated 
than it was assumed in BC, a fact which, as we showed in Sec- 
tion 3, seems to play an important role in the nature of transport 
in the phase space of Hamiltonian (1). 



5. Existence of irregular families 

In two previous papers (BC; Barbanis 1996) the authors investi- 
gated the bifurcation and the evolution of known periodic orbits 
belonging to some BL's. It was found that a small number of 
these families form pairs that do not have any direct or indirect 
connection with the periodic orbits of the unperturbed system. 
However, for a given multiplicity, there are many other bifur- 
cating families, which may play the role of connecting bridges 
between those pairs. 

Following these considerations, we studied the evolution of 
the bifurcations of some specific multiplicities from the basic 
orbits la and lb. We have selected the multiplicities 1 1 and 29 
for three reasons: (I) 1 1 and 29 are prime numbers, so that there 
are no families of these multiplicities resulting from lower mul- 
tiplicity families, except from those with multiplicity one. (ii) 
The number of bifurcating families is neither small nor great, 
(iii) There are a lot of orbits of these multiplicities on the vari- 
ous BL's and FTL's for e = 0.560 and 77 = 0.20. 

The detailed discussion in Appendix A and the examples in 
Appendix B lead to the conclusion that pairs of families with- 
out any direct or indirect connection to a basic orbit of the un- 



perturbed system (e.g. the six pairs appearing in Fig. A. 1 ) exist 



in a 3-D system. This happens mainly for high multiplicities, 
exactly as one would expect considering the analogous case of 
irregular families in 2-D systems. It should be noted also that, 
besides the above mentioned classical case of irregular fami- 
lies, in this work we have found a new kind of irregular fam- 
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Fig. 3. The contour lines of constant LCN at the value A = 
10~ 4 for a coarse mesh (Ax = Ay = 0.005, dashed line) and 
a fine mesh (Ax — Ay = 0.001, solid line). Regions of A > 
10~ 4 and A < 10 -4 are marked with a H and an L, respectively 

ilies, which we call a double pair, i.e a pair consisting of two 
single pairs (see Fig. |B.6] and its discussion in Appendix B). 

6. Summary and conclusions 

Barbanis and Contopoulos (1995), in studying the distribution 
of the p.c.'s of the periodic orbits in a well known 3-D model 
potential with the x — y plane, have discovered a noteworthy 
order. Except for a couple of gaps, the p.c.'s of periodic orbits, 
whose multiplicities form Farey-tree sequences, are arranged 
along continuous and smooth-looking lines. The p.c.'s along 
these lines have some fractal properties, arising from the fact 
that they belong to Farey-tree sequences of various orders. In 
the present work we tried to study the properties of these lines 
in association with the values of the LCN's of the area span by 
them. 

We have found numerical evidence that the contour lines of 
LCN's show fractal properties. We have also found numerical 
evidence that the splitting of initially continuous and smooth- 
looking BL's and FTL's, as well as the ensuing formation of 
gaps, is not an exception, as was implied in BC. As a rule, it 
is observed in regions of the x — y plane characterised by high 
values of the LCN's and appears not only in BL's, but in higher 
order FTL's as well. 

Since our results are only numerical and involve only a 
small number of periodic orbits, the above results cannot be 



considered as firm proof that the distribution of BL's and FTL's 
as well as the contour lines of LCN's have, beyond any doubt, 
fractal properties. However we feel that there is enough evi- 
dence that the above two geometrical objects (BL's and FTL's 
on the one hand and contour lines on the other) are closely 
related and that, furthermore, they show strong indications of 
fractal behaviour. Therefore we think that the study of some 
selected regions of the x — y should be done in more detail, 
by calculating LCN's in even finer mesh and by comparing the 
contour lines to higher order FTL's, in order to establish to a 
higher level of confidence the fractal nature of both geometri- 
cal objects. 

As far as irregular periodic orbits are concerned, we have 
found that the majority of the periodic orbits forming BL's and 
FTL's either bifurcate directly from la or lb or they are con- 
nected indirectly with them. However, by considering families 
of periodic orbits with large multiplicities, we find that some of 
them do not have any connection with la or lb but they form a 
pair with another family. In the present work we have encoun- 
tered also the more complex situation of a double pair, i.e. a 
pair whose members are pairs too. 

Acknowledgements. The authors would like to thank Prof. L. Martinet 
for a critical reading of the first version, which improved the final form 
of the paper. H.V. would also like to thank Prof. J.H. Seiradakis and 
Mr. K. Tsiganis for several helpful discussions. 

Appendix A: Search for irregular families; bifurcations 
from the orbit la 

Tables Al and A2 give the bifurcating families from the basic 
family la with multiplicities 29 and 1 1 respectively, when < 
e < 0.560 and r/ = 0.20. The corresponding bifurcations from 
lb are given in Appendix B. The bifurcations from la and lb 
take place when the index 



a = -2cos(2?ri?) 



(A.l) 



(where R is the rotation number of the bifurcating orbit) is 
equal to one of the two indices of the stability, s, of the or- 
bits la or lb (Contopoulos & Barbanis 1994). In this way we 
calculate ebif, i.e. the value of the parameter e (for -q = 0.20) 
for which one of the indices of the stability of the orbit 1 a or 
lb, calculated in quadruple precision, is equal to the index a. 

Each Table contains the rotation number, R, the index a, the 
values of ebif, the name of each bifurcating family and the way 
that this family evolves. The name of each bifurcating family 
from la or lb consists of a letter (a or b) and two numbers. 
The first number gives the multiplicity, to (11 or 29), and the 
second the nominator of R. For example, lla7 designates the 
family with multiplicity 11 which bifurcates from la and has 
R = ll\ 1 . Orbits of families whose nominator of R is odd have 
two p.c.'s with the (x — y) plane, while those with even nomina- 
tors have only one. In the case of one p.c. there are two different 
bifurcating families. We distinguish one of them from the other 
with an ac cent to its family name, i.e. 1 la6'. 

In Fig. AT we have drawn the BL's A and A', B and B', C 
and C'. The origin and the end of these lines are the orbits 4(+) 
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Fig.A.l. The lines A and A', B and B', C and C' are around 
orbit la. On these lines the p.c.'s of the periodic orbits with 
multiplicity 29 are depicted, when e = 0.560, r\ = 0.20. The 
orbits of the corresponding families of the bifurcations from 
orbit la of Table Al are also given, omitting from each name 
the prefix 29. 

and 1 1(-) for A and A', 3(-) and 2c(+) for B, 3'(-) and 2c for B', 
3(+) and 2c(+) for C and C'. For clarity reasons we omit the 
prefix 29 from the name of the families in Fig. 4.1 and Table 
Al. 

The known orbits of multiplicity 29 on the above lines are: 



: 29 and 29', 29A and 29 A' 
: (31, (32 and (32', /33, (34 and (34' 
cl,c2 and c2', c3, c4 and c4', c5, c6 



- On the lines A and A' 

- On the lines B and B' : 

- On the lines C and C : 
and c6' 

- On other FTL lines (not shown here) lie the termination 
points of the bifurcating families from la, with i?=5/29 to 
8/29, namely a5, a6 and a6', a7, a8 and a8' (see Table Al) 

The families of Table Al, which correspond to the first four 
rotation numbers, become highly unstable as e increases. This 
is the reason why the computations were not continued until 
e = 0.560. 

The corresponding orbits of the families with i?=9/29 to 
13/29, when e = 0.560, are the points cl, c2 and c2', c3, c4 
and c4', c5 of C and C' respectively. Similarly, the orbits of the 
families with i?=15/29 to 17/29 coincide with the points (31, (32 
and (32', (33. The families al4 and al4' pass through the point 
29b of the spiral Eb and 29a of Ea, respectively. The lines Ea, 
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0. 065 
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Fig. A.2. The family 29A forms a pair with the family 29 Aa. 
From the common point S (e im j n , 77) emanate the family 29A1 
and 29Aal by changing 77. The family 29 A 1 stops at the point 
SI (e im in,?7iX while 29 Aal atthepointS2 (e imin , 772)- No con- 
tinuation of 29 Aal has been found. On the other hand at the 
point SI two other families, i.e. 29 A2 and 29 A3 emanate with 
i] = 771 and changing e between e im i n < e < e 2m i n . The two 
families are connected at the point S3 (e 2m ; n , 771). No connec- 
tion of S3 with la has been found. 

Eb belong to the spiral structure around la as seen in Fig. [jj 
(see also Fig. 3 in Barbanis 1993). 

The families al8 and al8' are connected indirectly with the 
families (34 and (34' respectively thro ugh two other families. A 



similar example is shown in Fig. |A.3 . 

The family al9 bifurcates at eut = 0.085467 and exists for 
e < ebif until e m j n = 0.048765. At this value it is connected to 
a family which becomes highly unstable at e = 0.305. 

Each of the remaining orbits 29 and 29A of A , 29' and 
29 A' of A', c6 of C, c6' of C' belongs to a family which termi- 
nates at some minimum e, where a new family emanates, so that 
six pairs are formed. E.g. the orbit c6 is the 14th member of the 
basic sequence C. This family terminates at e m ; n = 0.53507. 
From this termination point emanates the family ca, w hich at 
e = 0.560 reaches the point ca on a FTL (see Fig. A. 1 ). 

Summarizing, the bifurcating families evolve, with respect 
to e, in the following way: 

a) The characteristics of many families move away from the 
original parent family and, for e = 0.560, they cross a BL 
or a FTL. 



8 B. Barbanis et al.: Imperfect fractal repellers and irregular families of periodic orbits in a 3-D model potential 

Table A.l. Bifurcations from la with m=29 



R 


a 


Cbif 


Fam. 


Comments 


1/29 


1.7JJ ill 


0467 1 4 


al 


fiicrhlv litiQt 


LI Ly 


1 s i ^ i ^osa 


n 1 (\aai(\ 


a2 
a2' 




Dt Ly 


1 ^QO 1 7 
- 1 . jyLLoOL J 


U, 1 J / jj / 






A/OQ 
^1 Ly 


-l.LyH 1 1 Lj 1 


u.zu /uzu 


a4 
a4' 




5/29 


-0.93681688 


0.253711 


a5 


jeetgtson 


6/29 


-0.53505668 


0.297548 


a6 
a6' 


> 


7/29 


-0.10827782 


0.338633 


a7 


> 


8/29 


0.32356399 


0.377069 


a8 
a8' 


> 


9/29 


0.740277631 


0.412949 


a9 


a9=cl 


10/29 


1.12237413 


0.446314 


alO 
alO' 


al0 = c2 
alO' =c2' 


11/29 


1.45199098 


0.477053 


all 


all=c3 


12/29 


1.71371435 


0.504489 


al2 
al2' 


al2=c4 
al2'=c4' 


13/29 


1.89530634 


0.525074 


al3 


al3=c5 


14/29 


1.98827591 


0.5184375 


al4 
al4' 


al4=29b 
al4'=29a 


15/29 


1.98827591 


0.467095 


al5 


al5=/31 


16/29 


1.89530634 


0.395884 


al6 
al6' 


al6=/32 
al6'=/32' 


17/29 


1.71371435 


0.310435 


al7 


al7=/33 


18/29 


1.45199098 


0.208203 


al8 
al8' 


al8^ 134 
al8'-> /34' 


19/29 


1.12237413 


0.085467 


al9 


al9^al91 



b) Some families terminate at a maximum e max or a minimum 
e m i n . In this case there are three possibilities: 
(I) Such a family is connected to another family bifurcat- 
ing also from the same orbit la or lb and terminating at 
the same e (see e.g. Ilb2' -^llb8', Table Bl). 
(ii) From the same parent orbit, 1 a or lb, bifurcate two fam- 
ilies, one of them terminating at e max and the other at 
£min- Another family starting at e m i„ terminates at e max 



(see Fig. A. 3). In few cases there are more than one in- 
tervening families, 
(iii) From the bifurcating family, which stops at e max , an- 
other family emanates which, for e = 0.560, crosses a 
BL or a FTL (see e.g. 1 lb3 reaches 1 If of F, F', Table 
Bl). 
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Fig.A.3. The two branches of the family lib connect to the 
branches of the family 1 lba, at the points a, b. At these points 
emanate the branches of the family b', which, by changing r\ 
from 0.20 to 0, reach the points al, bl. At these points termi- 
nate the branches of the family llab, which bifurcate from la 
whene = 0.1227395, T) = 

c) In bifurcating families with small rotation numbers the 
computations stop sooner than e = 0.560, because the or- 
bits of these families become extremely unstable. 

Some pairs of families that are connected at some e m j n have 
no direct connection with la. The question is: Is there some 
other family bifurcating from la, for some pair of values of e 
and r), which passes through the point from which emanates the 
pair when e = 0.560 and r\ — 0.20 ? 

Fig. A. 2 represents the pair of 29A and 29 Aa; their com- 
mon point, S, corresponds to e im i n = 0.5504835, r\ — 0.20. 
Keeping e = e im j n = const and varying n we find two other 
families, 29 A 1 and 29Aal, emanating from the common point 
S. The family 29 A 1 stops at the point, SI, corresponding to 
minimum 771 = 0.181232, while the family 29Aal stops at 
S2 when 772 = 0.184350. All efforts to find a continuation of 
29Aal were unsuccessful. On the contrary, from the point SI 
emanate the families 29 A2 and 29 A3 for 77 = rji = const and 
e varying from the value e im i„ to e 2m i n = 0.52840725, where 
the two families terminate at the common point S3. No connec- 
tion of S3 to la has been found. 

In the case of multiplicity 1 1 we find that all orbits connect, 
direc tly o r indirectly, to la, except for orbit 1 lb with two p.c.'s 
(Fig. A^ ). This orbit is a Farey tree orbit on the lines B and B' 
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Table A.2. Bifurcations from la with m=l 1 
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Fig.A.4. The p.c.'s of the periodic orbits with m=ll are illus- 
trated with the corresponding bifurcating families from la. All 
the orbits on the BL's are connected directly with la, except 
of 1 lb, forming a pair with 1 lba which is connected indirectly 
with la (see Fig A. 3) 
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Fig. A.5. The two branches of the bifurcating families Hal and 
lla7from lastopate max = 0.3099459 and e min = 0.1522655 
respectively. The family 1117 connects these families. 



R 


a 


Ebif 


Fam. 


Comments 


1/11 


-1.68250707 


0.138524 


Hal 


llal^lla7 


2/11 


-0.83083003 


0.265946 


lla2 
lla2' 


point lla2 
point lla2' 


3/11 


0.28462968 


0.373682 


lla3 


Ha3=llA of A, A' 


4/11 


1.30972147 


0.463427 


UaA 
lla4' 


Ila4=llc2 of C 
Ila4'=llc2' of C' 


5/11 


1.91898595 


0.527056 


lla5 


lla5=llcl of C, C 


6/11 


1.91898595 


0.409869 


lla6 
lla6' 


lla6=ll of B & A 
lla6'=ll' of B' 


7/11 


1.30972147 


0.155080 


lla7 


lla7->llal 



between the basic members 3(-) and 5(+). The two branches of 
family lib stop at e imin = 0.15204035, 77 = 0.20 at the two 
p.c.'s a, b (Fig. A3). From these points emanates the family 
llba which, for e = 0.560, passes through points llba (see 



Fig. A. 4). The two families form a pair with no direct connec- 
tion to la. However, there is an indirect connection to la. In 
fact, starting from the common points a and b with constant 
e = timin and diminishing 77 from 77 = 0.20 to 77 = 0, we reach 
points al andbl. Continuing this family, keeping 77 = and re- 
ducing now e until the value ebif = 0.1227395, we find that the 
new family llab bifurcates from la. Therefore, families lib 
and 1 lba do not form a pair of irregular families, because this 
pair is connected indirectly to la. 

Since our results are numerical, it should be useful to con- 
firm them through further work for the following two reasons: 

(I) The search for bifurcations from la and lb with multiplic- 
ity 29 has been confined within the interval < e < 0.560 
and 77 = 0.20 only. One may argue that, although highly 
improbable, bifurcations for e > 0.560 may play the role 
of connecting bridges of the pairs, considered here as irreg- 
ular families, to orbits la or lb. 

(ii) Following the above argument as well as for reasons of 
completeness, it would be desirable to find also the bifur- 
cations of la and lb for e = 0.560 and 77 > 0. 

Table A2 shows the bifurcations of the orbits with multi- 
plicity m=ll from la when < e < 0.560, 77 = 0.20. 

On the lines of Fig. A.4 we represent the known orbits with 
m=l 1, when e = 0.560, 77 = 0.20. These orbits are: 

- On A and A' : 1 1 (also on B), 1 1 A 

- On B and B' : 1 1 (also on A) and 1 1', 1 lb 

- OnCandC : llcl, llc2and llc2' 

- On other lines (not showing) 1 la2 and 1 la2', 1 lba 

The bifurcations of Table A2 are related to the above orbits 
as follows. 
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Table B.l. Bifurcations from lb with m=l 1 



R 


a 


Ebif 


Fam. 


Comments 


1/11 


-1.68250707 


0.168452 


llbl 


llbl=llbl 


2/11 


-0.83083003 


0.344563 


llb2 
llb2' 


llb2= lleofE 
llb2' -> llb8' 


3/11 


0.28462968 


0.525391 


llb3 


1 lb3=l If of F, F' 


8/11 


0.28462968 


0.290248 


llb8 
llb8' 


llb8— >lld of D 
llb8' -» llb2' 



Family Hal is connected to family lla7 through another 



family 1117, as shown in Fig. A. 5. Family Hal terminates at 



£max = 0.3099459. Family lla7 exists for e < ebif and ter- 
minates at e m i n = 0.1522655. Family 1117 starts at e m j n and 
terminates at e max . 

For e = 0.560, i] — 0.20 th e fo llowing families pass 



through the points illustrated in Fig. |A.4 



- Family 1 la2 passes though point 1 la2, lying on a FTL (not 
shown). Similarly lla2' passes through lla2' 

- Family 1 la3 passes through point 1 1 A on A and A'. 

- Family 1 la4 reaches point 1 lc2 on C, while 1 la4' reaches 
point 1 lc2' on C'. 

- Family 1 la5 reaches points 1 lcl on C and C' 

- Family 1 la6 passes through the cross point 1 1 of B and A, 
while 1 la6' reaches point 1 1' on B'. 

Families which start at points 1 lb on B and B' and 1 lba 
form a pair when e = 0.15204035, rj = 0.20. This pair has an 
indirect connection to la, as illustrated in Fig. A. 3. 



Appendix B: Search for irregular families; bifurcations 
from the orbit lb 

Table Bl gives the bifurcations with multiplicity m=ll from 
lb. Families with R = 4/11 to 7/11, 9/11 and 10/11 are not 
included, because the corresponding values of a are outside the 
interval of th e present study of < e < 0.560 and rj = 0.20. 

Fig. |jy] depicts the BL's E and E', F and F', D and D' 
and the places of the orbits with m=ll and m=29 when e = 
0.560,?7 = 0.20. 

The known orbits on these lines are: 

a) m=ll 

- OnEandE' : lie, lie' 

- OnFonF' : 1 If 

- OnDandD' : lid, lld'd 

b) m=29 

- On E and E' : 29el, 29el', 29e2, 29e3, 29e3' 

- On F and F' : 29f, 29f 

- On D and D' : 29dl, 29d2, 29d2' 



Fig. p. 2| depicts the characteristics of the families in Table 
Bl, except for family llb3, because its characteristic is very 
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Fig. B.l. The basic lines E and E', F and F' cross each other at 
the periodic orbit lb, while D and D' lie close to lb. On these 
lines we show some basic orbits with small multiplicities and 
the p.c.'s of the periodic orbits with m=l 1 and m=29. 



close to that of family llbl. Family llbl reaches points llbl 
when e = 0.560, while family 1 lb3 reaches the two points 1 If 
of F, F'. 

Family llb2 passes through point lie of E. Family llb2' 
reaches e m ax b8 , = 0.46578405, where it is connected to family 
1 lb8'. Family 1 lb8 reaches e maXb8 = 0.4622105, where fam- 



ily 1 1 8d emanates. This family connects llb8 to lid, which 
starts at e m i lld = 0.443468 and, when e = 0.560, it reaches 
point 1 Id lying on D. 

Family lid' terminates at e m ; n = 0.443836, where it joins 
family 1 le', which for e = 0.560 passes through point 1 le' of 
E' (Fig. B.2). Families lie' and lid' form a pair, not having 
any obvious connection to the basic orbit lb. 

Table B2 shows the bifurcating families with multiplicities 
m=29 from lb. As we mentioned before, the computation of 
orbits with small R, i.e. R = 1/29, 3/29 and 4/29, stops before 
e = 0.560 because these orbits become highly unstable. 



Fig. B.3 shows the indirect connection of family 29b2 to 
the family 29e3, starting from orbit 29e3 on E. Family 29b2 
terminates at e maX2 = 0.3339236, while family 29e3 at e m j n = 
0.2491375. The family connecting 29b2 and 29e3 is 292e. 

Family 29b2' is connected to 29b20' at e m ax 9 / = 
0.3334738. On the other hand family 29b20 stops at e = 0.360. 
Family 29e3', starting from point 29e3' on E', reaches e m ; n = 
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Fig. B.2. Location of the characteristics of the bifurcating fam- 
ilies of Table Bl parametrized by e (r\ = 0.20) until e = 0.560. 
Only the last points of family 1 lb3 are given. The pair of the 
families lie' and lid' which start through the corresponding 
orbits on E' and D' s also shown. Family 1 lb8 is connected to 
lid through the intervening family 1 18d. 



0.2491375, where family 29x emanates. The computation of 
29x stops at e = 0.408, r\ = 0.20. 



Fig. B.4 shows the connection of the branches of fami- 
lies 29b5 and 29b21 through two other families, i.e. 2951 and 
2921. Families 29b5 and 29b21 reach e maXB = 0.44357242 



and e n 



= 0.4545552 respectively. Families 2951 and 2921, 



which are connected at e m i„ = 0.4394852, start, respectively, 
in these points. 

For e = 0.560 family 29b6 reaches point 29d2 of the FTL 
sequence 5(+), ....29d2(-), 24(+), 19(0, 14(+), 9(Q on D, while 



29b6' is going through point 29d2' on D' (Fig. B.l). Similarly 
for e = 0.560 family 29b7 reaches points 29dl of the basic 
sequence D and D', namely the sequence with first member the 
orbit 5(+) and increment 4. 
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Fig. B.3. The bifurcating family 29b2' is connected to 29b20'. 
Families 29b2 and 29e3, corresponding to point 29e3 on E, are 
connected indirectly through family 292e. Family 29b20 be- 
comes highly unstable at e = 0.360. Family 29e3' forms a pair 
with 29x. 

Families 29el and 29f, going through points 29el on E 
and 29f on F respectively, form a pair connecting at e m i n = 
0.5520015. 

The two branches of family 29e2, which start from points 
29e2 on E and E' when e = 0.560, parametrized by e (77 = 
0.20) stop at point SI when e millo2 = 0.43790 (Fig. ^|). A 
new family 29e21 emanates here and terminates at point S 
when e max = 0.4437953. This pair is connected to another pair 
which is formed by two families having two branches also. The 
first family of this pair starts at e = 0.560 and -q — 0.20 from 
points 29g, lying on FTL's and terminates at point S2, when 
emin.g = 0.4435759. At this point emanates the second fam- 
ily 29gl of this pair. This family terminates at S together with 
29e21. Such a double pair, having no connection to an orbit of 
the unperturbed system, is noticed for the first time. 
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Fig. B.5. The direct connection of 29b8 to 29b22 as well as the 
indirect connections of 29b22' to 29f' through family 2922f' 
and of 29b8' to 29el' through family 298el' are illustrated. 
Families 29el and 29f, starting from points 29el on E and 29f 
on F respectively, form a pair which has no obvious connection 
to lb. 



B. Barbanis et al.: Imperfect fractal repellers and i 

t 1 1 1 1 1 1 1 




-0.10 - 



0.04 0.08 x 0.12 

Fig. B.6. Two pairs of families are connected at the common 
point S. The first pair, consisting of families 29e2 and 29e21, 
has as common point S 1 . The second pair of the families 29g 
and 29gl is connected at S2. 
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Table B.2. Bifurcations from lb with m=29 



R 


a 


Cbif 


Fam. 


Comments 


1/29 


-1.95324111 


0.063890 


29b 1 


highly unst. 


2/29 


-1.81515084 


0.127362 


29b2 
29b2' 


29b2^29e3 
29b2' ^29b20' 


3/29 


-1.59218613 


0.192211 


29b3 


highly unst. 


4/29 


-1.29477257 


0.258516 


29b4 
29b4' 


> 


5/29 


-0.93681688 


0.326001 


29b5 


29b5^29b21 


6/29 


-0.53505668 


0.394292 


29b6 
29b6' 


29b6 = 29d2 
29b6' = 29d2' 


7/29 


-0.10827782 


0.462973 


29b7 


29b7 = 29dl 


8/29 


0.32356399 


0.531620 


29b8 
29b8' 


29b8^29b22 
29b8' ->29el' 


9/29 


0.74027631 


e b if>0.56 






20/29 


0.74027631 


0.078653 


29b20 
29b20' 


29b20 stops 
29b20' ^29b2' 


21/29 


0.32356399 


0.270218 


29b21 


29b21^29b5 


22/29 


-0.10827782 


0.516148 


29b22 
29b22' 


29b22^29b28 
29b22' ->29f 



